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Measurements of the r\ meson production with a polarised proton beam in the reaction 
pp — > ppr\ have been carried out at an excess energy of Q = 40 MeV. The dependence of the 
analysing power A y on the polar angle 9* of the n meson in the center of mass system (CMS) 
has been studied. The data indicate the possibility of an influence of p- and d-waves to the close to 
threshold r\ production. 
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I. INTRODUCTION 

Several measurements on the r\ meson production in 
the proton-proton interaction covering a 100 MeV excess 
energy range were performed at different accelerators. 
The determined total cross sections |], 0, j| 0, || |) , as 
well as their differential distributions Pf pf H, TO] trig- 
gered intensive theoretical investigations aiming to un- 
derstand the production mechanism on the hadronic and 
quark-gluon level. 

In the theoretical descriptions of the ^-production in 
nucleon-nucleon collisions the excitation of the Sii(1535) 
resonance plays a decisive role. The hitherto performed 
studies with the aim to describe the total cross section 
show a dominance of this virtual S\\ nucleon isobar in 
the close-to-threshold production of the r\ meson. The 
excitation of this intermediate state results from a one 
meson exchange (e.g. n, rj or p) between the two nucleons 
followed by a strong coupling of the 77.ZV system to the S\\ . 

Near threshold the energy dependence of the total cross 
section results from a three-body phase space modified 
by a strong nucleon-nucleon final state interaction and 
a significant contribution of the attractive interaction in 
the rjp system. Since several existing models with differ- 
ent scenarios of the excitation describe the existing data 
well, a confrontation of the predicitions with other ob- 
servables is needed in order to distinguish between them. 
The measurements with polarised beam should settle the 
on-going discussion whether the rj production is domi- 
nated by p 0, |R| [lij, w |Jj or V El exchange. 
The interference between considered amplitudes causes 
a different behaviour - depending on the assumed sce- 
nario - e.g. of the r\ meson angular distributions. These 
differences are too weak in the close-to-threshold region 
to discriminate between different models. Yet, the pre- 



dictions of the analysing power depend crucially on the 
assumed mechanism |l7], [l8j . 

So far the only measurement of the analysing power has 
been performed Q at an excess energy Q — 1805 MeV. 
In the present experiment, the analysing power close to 
the production threshold is determined and results for 
the interference terms from contributing partial waves 
are presented. A comparison with theoretical predictions 
will be discussed in section |v|. 

Section || contains the description of the experiment 
and the method to extract the pp — ► ppr\ events. The 
following section introduces definitions and gives a theo- 
retical overview. In section [V, the results are presented. 



II. EXPERIMENT 

Measurements of the pp — > ppr\ reaction were per- 
formed at the internal experiment COSY-11 at the 
COoler SYnchrotron COSY |20) in Jiilich with a beam 
momentum of pbeam = 2.096 GcV/c corresponding to an 
excess energy of Q = 40 MeV. During the experiment, 
cycles of about ten minutes for the two different beam 
polarisations were adjusted. 

Using a hydrogen cluster target [|l| in front of one of 
the regular COSY dipole magnets, the experimental facil- 
ity acts like a magnetic spectrometer. Positively charged 
particles in the exit channel are bent towards the inte- 
rior of the ring where they are detected in a set of two 
drift chambers. Tracing back the reconstructed trajecto- 
ries through the magnetic field to the interaction point 
allows for momentum determination. Particle identifi- 
cation is achieved by a time of flight measurement over 
a distance of 9.4 m between two scintillation hodoscopes. 
For further details, the reader is referred to reference [191. 
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Figure 1 (a) shows that the method allows for a clear sepa- 
ration between pions and protons and hence for the iden- 
tification of events with two protons in the exit channel. 
The not registered particle system X - either a sin- 




(a) (b) 

FIG. 1: (a) Invariant mass spectrum for events with two re- 
constructed tracks. Besides the clear proton peak a second 
signal stemming from pions is observed, (b) Missing mass 
squared for events with two protons in the exit channel. Lit- 
erature values J5^] for particle masses are indicated by arrows. 

gle meson (the r\ in the present case) or a multi me- 
son system - is identified by means of calculating its 
mass m\ = (P beQm + P tar get -Pi - P2) 2 , while F beam 
and Ptarget denote the four momentum of the beam and 
target proton in the initial channel and Pi, P2 those of 
the two registered protons. The missing mass spectrum 
for events with two identified protons is shown in Figure 



1(b) for the entire beam time. Besides the clear 77-signal 
there is obviously a 7r°-peak resulting from the reaction 
pp —>■ ppir°. Furthermore, a broad yield due to multi pion 
events with the lower limit given by m x = (2m 7V ) 2 and 
the upper limit by m 2 x = {y/s- 2m p ) 2 = 0.345 GeV 2 /c 4 
is observed. The increasing event rate towards higher 
missing masses is due to the higher acceptance of the 
COSY-11 detector for two protons with small momenta 
in the center of mass system (CMS). 

The monitoring of the geometrical dimensions of the 
synchrotron beam and its position relative to the tar- 
get [^3) enable to achieve a mass resolution of cr mj) = 
1.6MeV/c 2 . The much broader peak of the ir° is due to 
the error propagation which worsens the mass resolution 
with increasing excess energy || . 



III. GENERAL DESCRIPTION 
A. Definitions 

A detailed theoretical derivation of the analysing power 
was recently published for the case of the pp —> ppir re- 
action g5| |2(| . For the 77 production the description is 
analogue since in both measurements the initial channel 
is fixed to isospin 1=1. Therefore, the different quantum 
numbers for 7r° (as a member of an isotriplet) and the 
isoscalar rj are irrelevant. 

In the given experimental situation a convenient choice 



of the three axis is: 

Pbeam 



111 



'beam \ 



P A y x z 
-=r and x = 7— -r, (1) 

|P| \yy-A 



where P indicates the polarisation of the COSY beam. 

In the COSY-11 experiment, the two four momenta 
of the final protons Pj = (E*,p*) are measured. The 
CMS momentum of the rj meson is q = —(pi +P%)- The 
proton momentum in the pp rest-system is denoted by p 
For later purposes, Figure |2| depicts the definition of the 
used polar- (9) and azimuthal angle (ip). The indices p 
and q will refer to the pp rest-system and the rj meson 
in the CMS, respectively. The angle 9 P will be choosen 
such, that < p < ir/2. This choice guarantees that all 
observables are invariant under the transformation p — > 
— p as required by the identity of the two protons in the 
final state. 





X 




IV 




-->tP 


1 1 / 




1 






, — / 


NcV 




/ 

/ 


/ / ' 




/ 


// 




/ 


fir 


/ \ 


/ 



beam 



e 



FIG. 2: Definition of the angles. 9 is defined as the angle 
between momentum vector and the z-axis, ip between the x- 
axis and the projection of p onto the x-y-plane. 



B. Observables 

The differential cross section for a reaction with a po- 
larised beam is given in terms of Cartesian polarisation 
observables by p7fl 



*(0 = *o(0(i + X>i-A(0 



(2) 



i=l 



where Pi and Ai denote the beam polarisation and the 
analysing power in the given reference frame, oo(£) indi- 
cates the total cross section in case of no polarisation. In 
the upper formula, the abbreviation 



*(0 = 



d 3 <T 



(0 



is used where £ denotes the set of the five variables which 
are kinematically completely describing the exit channel, 
namely (9 p , (p p , 6*, <p*, E pp ). The kinetic energy E pp 
of the two final protons in their CM system is given by 

with y/su — 2y Ip 1 + m p as the energy 



'S12 — "i»p 
in the pp subsystem 
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In the given case of the general experimental condi- 
tions, the beam polarisation is - due to the magnetic 
fields in the accelerator - forced to be P = (0, P y ,0) T 
and hence formula (|2|) simplifies to 



(3) 



The asymmetry e - obtained from the difference in the 
yields with beam polarisation up and down - defined by 



Nt-N 



7V T + N t 



(4) 



forms the basis for deducing the analysing power while 
Nf (N±) denote the experimental number of events for 
spin up (down). With known luminosity C, efficiency £ 



and measured time dt, dN is related to the cross section 
by dN^i = £ ■ £ ■ a] _i ■ dt. In combination with equation 
(0), one can deduce from (Eh that 



MO 



C-rel • N-t 



AT, 



N r P r - C rel ■ 1\ ■ A 



(5) 



where the relative time-integrated luminosity C re i is de- 



fined by C re i ■— 



fCydtj. 



In equation M, the efficiency 



cancels out because of the independence on the spin as 
long as the bin size of A£ is small enough so that the 
efficiency can be assumed to be constant. 



With the definitions given in section III A , the angular 



dependence of the spin-dependent cross section can be 
written as 



<ro(OMO = {[Gf + Gf (3 cos 2 1)] sine; + [Hf + P° + Hf {3 cos 2 9 p ~ 1)] sin 2#* } cos <^ 
+ [Hf + K y0 + Gf cos 9* q + iff °(3 cos 2 0* - 1)] sin 29 p cos ip p 

+ (Gf sin 9* + Hf sin 29* q ) sin 2 9 p cos(2^ p - ip*) + Hf sin 29 p sin 2 9* cos(2ip* - ip p ) (6) 

I 



The appearing literals 1 denote interferences of partial 
wave amplitudes. The relative angular momentum of the 
two outgoing protons in their rest system is denoted by 
capital letters l p — S,P,D..., the one of the r\ meson 
in the CMS by small letters l q = s,p,d..., while the 
usual spectroscopic notation is used. With this defini- 
tion, the single terms G^°,Hf,I y0 and K y0 correspond 
to (PsPp), (Pp) 2 , (SsSd) and (SsDs). 



The relative luminosity C re i — ^, as • T ,,, aa was ex- 

tracted via the elastic proton-proton scattering. To de- 
termine the elastic cross u^ s section according to equa- 
tion (pj) the analysing power was taken from ^9|. With 
the number of events N^ a resulting from the elastic pp- 
scattering C re i was calculated according to the definition 
given above: 



IV. RESULTS 

In order to extract the assymetry from the measured 
spinup and spindown events one needs the relative lumi- 
nosity and the average beam polarisation. 

Via a simultaneous measurement of the proton- 
proton elastic scattering at the internal experiment 
EDDA [^H, |29) the polarisation was determined for two 
time blocks 2 : 





time block 1 


time block 2 


P T 


0.381 ±0.007 


0.497 ±0.006 


Pi 


-0.498 ±0.007 


-0.572 ±0.007 



1 The superscript yO indicates a beam polarisation along the y-axis 
and an unpolarised target. 

2 The significant increase of the polarisation from the first to the 
second block is caused by improved tuning of the beam with 
respect to polarisation. 





time blockl 


time block 2 


Crel 


1.004 ± 0.004 ±° °° 2 


0.949 ±0.004+° °°! 



An integration of equation (g) over cos 9 p and ip p leads 
to the disappearance of several terms provided the exper- 
imental angular distribution covers either the full phase 
space with a constant detector efficiency or symmetri- 
cal ranges. Figure |^ shows the angular distributions of 
pp?7-events from a Monte-Carlo simulation which are nei- 
ther symmetric around 90° in case of tp p nor constant for 
both angles cos9 p and (p p . Therefore, the evaluation of 
the analysing power requires an efficiency correction. To 
correct the data the efficiency £(cos 9 P , (p p ) is determined 
via Monte-Carlo simulations. Using a GEANT-3 code for 
each event a detection system response was calculated 
and the simulated data sample was analysed with the 
same programme which is used for the analysis of the ex- 
perimental data. Weights w(cos9 p , <p p ) = l/£(cos9 p ,Lp p ) 
were applied during the final analysis of the experimental 
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data and hence the corrected number of events reads: 



i k 

E E w hk 

i k 



(7) 



while i and k run over the bins tp p and cos 9 P , respectively. 
The error is deduced with AAT^fc = y/N ijk to be 




AN cor = 



whereas the error of Wi t k was neglected because of a much 
higher statistic for the Monte Carlo simulations, so that 



An 



The influence of the strong proton- 

m,k N i]k o t 

proton final state- interaction (FSI) was included via the 
description with a Jost-function fl30j . Former acceptance 
studies on the dependence on the various Jost function 
prescriptions showed a change of the result of maximum 
10% |3lj, ||^]. An extensive discussion on the influence of 
the FSI reflecting itself in the density distribution of the 
Dalitz plot is given in |33|. Concerning that the FSI is 
known up to an accuracy of around 30% one can conclude 
that the upper limit for the total contribution to the er- 
ror is approximately 3% which is negligible compared to 
the high overall error of this first data sample. 
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FIG. 3: Angular distribution for the two proton angles cos# p 
and tfip obtained via Monte-Carlo simulations. 

Only after an efficiency correction one can remove 
the dependency from proton-coordinates in the analy- 
sis which is then the same as an integration over these 
variables so that equation (^|) simplifies to: 



d 2 a 



dilpdQq 



(0 A y (£)d cos 6p d(p p 



2tt Gf sin 9* + (Hf + F°) sin 29* cos 



(8) 



Due to the restricting dipole gap ip* is dominantly peaked 
around 0° - quite similar to the (p p distribution - but 
with a negligible peak around 180° which is not shown 
here but verified with MC simulations. Therefore, the 
analysis was performed with one single </?*-bin around 
±30°. Hence, equation (||) leads further to the separa- 
tion of the (PpPs)-interference (G\°) and the (Pp) 2 - and 



(SsSd)-terms (Hf and I y0 ): 



G\ 



i/0 



/(cos#*) dcos9* , 



(9) 



H f + I yo = 2_ J f (cos e* ) cos 0*d cos 9* , 
with/(cos0*)= / JfJ^r^AyiOdcosOpdtppd^ 



Defining 3 A^(cos6»*) := /// iV cor - (0 dQ p dtp*,, it is 
straightforward to show analogue to equations (|j) and 
(p|) that the integrated analysing power defined by 



A y (cos9*) :=/(cos0*)/ 



da 



d cos 9* 



can be determined via 
A y (cos9*)= £ " 



■N, 



n,1 



N, 



'i-l 



Pi ■ u 



1,1 



- (cos ( 



(10) 



(11) 



The calculation of A y needs the determination of the 
absolute pp — > ppr] events N v ^± in dependence of cos6>*. 
In section ||, the selection of the ppr) was discussed. The 
analysis was performed with 4 bins in cosO* starting at 
cos9* = —0.75 with A cos 9* — 0.5. A re presen tative 
missing mass spectrum is shown in Figure 4(a) where 
the background is fitted by a polynomial function. From 
this spectrum the number of events Nri+b including back- 
ground and 77-event are extracted. Subsequently, this 
background is subtract ed an d the number of events N v 
are determined (Figure 4(b) ). 
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FIG. 4: Spectra of the squared missing mass for events with 
two identified protons (a) with and (b) after subtraction of the 
background. The number of events was corrected according 
to equation (0). 

For the two time blocks, an error weighted mean value 
for A y is calculated. Figure |5| shows the analysing power 



3 In the following, limits of the integrations will be omitted as they 
are always the same. 
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i/O 



as a function of cos 8*. The extraction of G\ 

Hf + P° with equations (§) needs according to 
the knowledge of 

fact that tp q G [— g-, -g-j ™iuuw «^ ^- cos 

/„ 6 dip* which is due to the isotropy of the 
section in ip 



l COS 

7T 7T1 

' 6 ' 6 J 



which was taken from |7J. 
was considered with 



and 

© 

The 



da 



dcos 9* 



7r da 
3dn* 



cross 



(12) 
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FIG. 5: Dependence of the analysing power on the center of 
mass polar angle ff" q of the 77 meson. 



V. COMPARISON WITH THEORY 

The present data on the 77 meson production in 
nucleon-nucleon collisions referred to in section | show 
not only the 3-body phase space Q 2 -dependency and a 
modification due to the nucleon-nucleon final state in- 
teraction but also a significant influence of the nucleon- 
meson interaction in the case of the r\p system. As men- 
tioned above, several models describe the existing data 
quite well although they are based on different assump- 
tions for the excitation mechanism of the Su resonance. 
For instance, Batinic et al. [jl6) or Nakayama, Speth and 
Lee Jl8| found a dominance of ir and 77-exchan ge i n the 
analysis of pp — ► pprj while Faldt and Wilkin |17|] con- 
clude a dominant p-exchange. 

Polarisation observables may be the right tool to dis- 
tinguish between the different models. Calculations for 
the analysing power in the reaction pp — > pprj show dif- 
ferent results depending on the underlying assumption 
for the one meson exchange model. Figure H presents 
results taken from references JTt) (dotted line) and 
(solid and dashed lines) for Q = lOMeV and 37MeV. 
The authors of the latter reference conclude in the full 
model calculations a dominance of 7r and ^-exchange 
(solid line). The dashed curve represents a vector dom- 
inance model with an exclusion of tt and ^-exchange for 
exciting the Su resonance. The triangles are the experi- 
mental results. 



The averaged values of A y and the cross section used 
for the integrations in equations (JsJ) are presented in ta- 
ble § 



COS q 


Ay 


w d "n [Mb] 

a COS 6 q J 


-0.75 ±0.25 


0.19 ±0.21 


0.31 ±0.01 


-0.25 ±0.25 


-0.02 ±0.09 


0.50 ±0.01 


0.25 ±0.25 


0.05 ±0.06 


0.50 ±0.01 


0.75 ±0.25 


-0.05 ±0.06 


0.31 ±0.01 



TABLE I: Analysing power as a function of the emission angle 
8 q of the r} meson in the CMS and the differential cross section 
obtained from Q with equation (p"2|). 

Finally, the integrations of these values 

1 tt 2 ^ dcosO* u q 
cos e* q 1 

Hf ± F° = \ y -j^-Ay cos 0; ■ A cos 6* . 
1 7T 2 ^ dcos9* u q q 



cos e* ? 



result in 



and 



Gf = (0.003 ± 0.004) 11b 



Hf + I y0 = (-0.005 ± 0.005) 11b . 
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FIG. 6: Analysing power for the reaction pp — » pprj in de- 
pendence on 8* for the two excess energies Q = 10 MeV and 
37 MeV. 



The observable structure of the experimental values 
show a slight deviation from the sin 6 q cos # g -dependence 
of both models, ft seems that the data favours the vector 
dominance exchange models. The more or less strong 
difference in the angular dependency of A y results from 
a vanishing Gf in both references. As this corresponds 
to the (PpPs)-term, an influence of the P-wave must be 
suspected but right now the experimental result for Gf is 
compatible with zero. A non-zero Gf would imply that 
Hf - describing the (Pp) 2 interference - should have 
a non negligible contribution, too. For further detailed 
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studies the data are not yet precise enough to disentangle 
the sum of Hf° and I y0 . At this time the results indicate 
the possibility of an influence of p- and d-waves to the 
close to threshold r\ production. 

VI. CONCLUSION 

The reaction pp — > ppr] has been studied at an excess 
energy of Q — 40 MeV. The final state has been kinemat- 
ically completely reconstructed and the analysing power 
has been determined. Qualitatively, the data seem to 
favour the calculations with dominant vector meson ex- 
change but definitive conclusions cannot be drawn due 
to the large uncertainties of the data. To allow a more 
rigorous comparison with theoretical calculations higher 
statistics experiments are required and already scheduled 
for 2002 at COSY-11. 
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